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. , Abstract 

> . 

We present a new method for reducing the Fredholm determinant associated with an underlying 
Birman-Schwinger operator to a finite dimensional determinant. Moreover, we compute explicitly the 
connection between the Fredholm determinant and the Evans function for travelling wave problems of 
3J, ■ all orders in one dimension. 

^j '. 1 Introduction 

The purpose of this paper, is to determine the spectrum of a linear differential operator L = Lq + v(x), 
where Lq is a constant coefficient differential operator of order n > 2, and v(x) is the Jacobian evaluated at 
the travelling wave <f>. The approaches adopted for solving such a problem in this paper are essentially, the 
Fredholm determinant and the Evans function. Both of these are analytic functions whose zeros coincide in 
location and multiplicity with the eigenvalues of the operator L. The Evans function or the miss-distance 
function [35], as introduced in [5] and extended for instance in [TJ, is defined as the Wronskian of the 
set of solutions u^ of the associated problem (L — X)u — which decay at x = ±oo, respectively. On the 
t-H , other hand, the Fredholm determinant is the determinant of a trace class operators that differ by the identity 

operator. Among its applications, we cite for example from mathematical physics [331 HU 123 [Ml QUI E3 ■ The 
motivation of this article is based on the fact that, numerically we have observed that the eigenvalue A = 
when it is embedded in the essential spectrum is a pole for the Fredholm determinant in some situations. 
Consequently, this implies that the Fredholm determinant is a rational function. So the behaviour of the 
Fredholm determinant near A = depends entirely on its rational expression. It turns out that under some 
• • . assumptions, the analytic rational expression is nothing but the Evans function and the Wronskian of the 

set of solutions corresponding to (Lq — \)u = 0. Hence the Fredholm determinant, which is an infinite 
determinant, is equal up to a non-zero analytic function of the spectral parameter, to a finite dimensional 
determinant. This observation is not new, since it has been proven in [33J [T21 H3J [T71 [HJ [2Z1 [Ml HI]: where 
the associated integral operator is either trace class or Hilbert-Schmidt operators. However, our approach of 
achieving similar result is different, as it uses the interpretation of the Evans function that is, the measure of 
the linear dependence of subspaces E^ defined by the set of solutions of (L — \)u = decaying at ±oo. By 
this, we mean that if the subspace E~ decaying at — oo (respectively +oo) is orthogonal to a given subspace 
decaying at +oo (respectively — oo) defined by the set of solutions of the adjoint problem, then necessary 
we must have that E~ and E + are linearly dependent. By doing so, we obtain a Gramian matrix D(\) 
of dimension less than n, which we call the matrix transmission coefficient by analogy to the transmission 
coefficient in Scattering Theory. The basic details of our reduction are as follows. 



m 



C3 



1. Under some assumptions, we prove that the determinant of D(\) is equal to the Fredholm determinant 
associated to a trace class Birman-Schwinger operator. 
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2. We observe that the Evans function E(X) is written as the product of the determinant of the matrix 
transmission coefficient and the determinant of the set of solutions of (Lq — A) = 0. This enables us to 
relate the determinant of the variable coefficient differential operator L — A to the Evans function. Thus 
we can define the Fredholm determinant as the quotient of the determinant of differential operators 
which are Fredholm. This implies that we can extend the notion of the Fredholm determinant to 
non-compact operators. 

3. Given the decomposition of the Evans function in [21 we extend the connection between the determinant 
of D(\) and the modified Fredholm determinant associated to a general first order system Birman- 
Schwinger operator [T5J |H] . 

The organisation of the paper is as follows. In Section 2 we construct the Green's function for our problem. 
In Section 3 we introduce the Birman-Schwinger operator corresponding to our model, and some of its 
properties. In Section 4, we make the connection of the Fredholm determinant to the determinant of 
-D(A). This connection enables us to extend the connection to integral operator belonging to weaker class 
than that of trace class. In section 5, we construct the appropriate Fredholm determinant for the case 
linia^-i-oo 4>{x) = <f> with (f> + ^ <p~ . This is done by constructing an equivalent Fredholm determinant such 
that its associated integral operator contains only the essential information of the original problem. Finally 
in Section 6, we provide some concluding remarks and outline future directions of investigation. 



2 Green's Functions 

A differential operator A defined on L 2 (R, C) with domain dom(„4) = L 2 (R, C) is invertible if it is injective 
and if its range, i.e Im(.A), is the entire space L 2 (R, C). This because we wish A to have both left and right 
inverses with these being the same. Assume A is invertible, then we write 



A~ x u(x) = \ g{x,y)v(y)dy, Vt; e L 2 (M,C), 



where g(x,y) is the Green's function. Throughout this paper, by invertibility we mean the existence of the 
left and the right inverse, i.e 

A~ l A = AA- 1 =I L 2. 

Let C be an nth order differential given by 



d" d 

C = a n (ir)-—- + . . . + oi(»)— + Oo(x). 

Theorem 2.1. f28f The Green's function g{x,y) satisfies Cg{x,y) = 5(x — y) in the sense of distributions 
if and only if : 

1- £g{x,y) = 0forx^ y; 

2. g(x,y) is n — 2 times continuously differ entiable in x at x = y; 

-\ f * x=y 



3. 



dx" 



— ^r (jump condition) . 



If the Green's function exists, then the solution of Cu(x) — f{x) is given by 



a 



b 



(x) = Kf(x) = / g(x,y)f(y)dy 



a 



where K denotes the inverse of the operator C. 



The maximal operator £ max corresponding to C is such that 

dom(£ max ) = {ue AC (I) n L 2 (I) : Cu G L 2 (I)}, 

where I CM. and AC(I) denotes the set of absolutely continuous functions on /. 

Proposition 2.2. \17\ Chap XIV. 3] The maximal operator corresponding to C and any interval is closed. 

The Green's function of a closed ordinary differential operator on a domain / is semi-separable, this means 
that g{x,y) has the following form [17, Chap XIV. 3] 

\Fi(x)Gi(y), y<x, x,y G I; 
[F 2 {x)G2{y), x < y, 

where F 3 {x) G L 2 (I,C nxn i) and G 3 {x) G L 2 (7,C^ xn ) for j = 1,2 and some m,n 2 G N. 
We recall a bounded operator £ : X — !> F acting between Banach spaces X and y, is Fredholm if its range 
is closed and the dimension of its kernel and cokernel are finite. Consequently we have, if K : X — > X is a 
compact operator then the operator I + K is Fredholm, where / is the identity map of X to itself [T7] . 
Consider the differential operator L obtained from linearising a nonlinear partial differential equation (PDE) 
on K about its travelling wave solution <fi(x). The operator L is defined on L 2 (R,C) and it is given by 

L = L Q +v(x), i£l, (2) 

where 

d™ d™ -1 d 

Lo= d^ +a "" l(a;) d^ rT + --- + ai(x) d^ +ao(x) ' XGM ' (3) 

and v{x) = V((f>(x)) where V is the Jacobian corresponding to the nonlinear part associated to the PDE. 

In what follows, we assume that (fy- l > is in L 1 and satisfies also 

lim *<«>(*) = (J"' l = °; (4) 

x|->oc I (J, ( > 1, 

where <p^{x) — <A l /dx l for some I < oo, and <j}°° G R. For simplicity, we take 0°° = 0, however the analysis 
that we will carry out remains the same for </>°° ^ or more generally when 0~°° ^ </> + °°, under the 
assumptions described in this paper. In that case, (Lq + V{<f>°°) — XI l 2 ) an d iy{x) — V(<fi°°)) substitute for 
(Lq — XI 1,2) and v{x), respectively. The eigenvalue problem reads 

{Lu = Xu 
w(±oc) = 0. ( ' 

Our objective is to determine the values of A for which the operator (L — XI^) : _ff"(R,C) — ► L 2 (R, C) is 
Fredholm of index zero. To this end, let's make the following hypothesis. 

Hypothesis 2.1. Assume that the aj are constant for i = 0, . . . ,n — 1, and the characteristic polynomial 
V n associated to the operator Lq — XI L i 

V n (\x) = / u" + a„-i^™- 1 + ... + ai/i + a - A = 0, VA G Q, (6) 

has k roots nf = /j,f(X), i = 1, . . . , k with Re(ft^~) > and n — k roots k~[ = /i^(A), i = k + 1, . . . ,n with 
Re(K~) < for a suitable Q C C which may contain the right-half complex plane. 

Under Hypothesis 12. 1[ the operator Lq — XI L 2 is Fredholm, therefore it has a bounded resolvent for A G O. 
Let K{X) : L 2 (R,C) — > L 2 (R,C) denotes the resolvent operator corresponding to Lq — A/^2, that is 

(K(X)u) (x, A) = / g(x, £, A)u(fc A) d£, Vu e L 2 (R, C). 
Jr 

Thus from Proposition 12.21 Theorem 12.11 and (JXJ) , we have the following definition. 



Definition 2.3. For A 6 il, we define the Green's function g(x, £, A) £ L 2 (R 2 , C) of the operator Lq — XI L 
by 



^a l (A)e K ^ a; -«\ x<£, x,(el; 

£ ai(\)e<(*-t\ £<x, 
where c\ii(X), i — 1, . . . , n are solutions of the following matrix equation 



(7) 
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l k+l 



-1 
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+\n— 1 



(4) 



+ \n-l 



< K k+l) 



n-1 



-K) n -V 



a 2 (A) 



WwJ 



f°\ 



\-y 



(8) 



Let detf denote the determinant of a Fredholm operator [23 HI] - Under our assumptions, K(X)v(x) is 
compact. Then by standard argument invoking the Fredholm alternative. We have A is an eigenvalue of the 
operator L if only if 

det F (7 L 2 + K(X)v{x)) = 0. (9) 

Thus A is an eigenvalue of L if the determinant of the Fredholm operator / + K{X)v{x) is equal to zero. In 
other words if —1 is an eigenvalue of K(X)v{x). 

Rewrite the eigenvalue problem ([5]) as a first order system of ordinary differential equation, that is 

Y' = (A(X) + R(x)) Y (10) 



where Y = (u(x), . 



,.("-!)! 



A(X) = 
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— On-1/ 



(11) 



and 



for some fixed m £ {0, . 



R( X ) = Z' 0(n-l) X (n-l) ((n _ 1)xl) \ (12) 

\VC<PWJ(lx(m+l)) 0(lx(n-m-l))/ 

, n — 2}. Let the constant-coefficient differential operator La be defined by 



La :H 1 (R,C n ) -> L 2 (K,C") 

y i-> dr/dx - A(A)r 

Since the Fredholm properties, and the Fredholm indices of Lq — XI^ and d/dx — A(A) are the same [U|, then 
the constant-coefficient operator La is Fredholm as well, for all AgO. In other words, the Fredholmness of 
the operator La is a consequence of the existence of an exponential dichotomy on R (see. [2]). We write 
K (X) : L 2 (R, C") -> L 2 (M, C"), the inverse of the operator L^ given by 

(Ko(X)Y) (x, X)=J G(x, C A)Y (£, A) d£, V y e L 2 (R, C"), 

where G(x, £, A) G L 2 (M 2 , C™ xn ) is the Green's matrix satisfying 

L A G(a;, £, A) = 0, for x ^ £, (13a) 

G(e + ,e,A)-G(r,c,A) = /„, (i3b) 

G(x, x_, A) — G(x,x + , A) = I n , (13c) 

G(±oo,£,A) =0. (13d) 

Using equation (|T3"|) . we define the Green's matrix G(x,£, A). 



Definition 2.4. For A G 0, we define the Green's matrix g(x, £, A) G L (R , C rixrl ) o/ </ie operator La by 



G(x,£,A) 



-r -(o:,A)z + (e,A), x<e, z,£g 

r + (x,A)Z -(C,A), £<z, 



(14) 



w/tere Y ~(x,A) eI 2 (l,C xi ), Y+{x,X) g i 2 (R,C Tlx ("- fc )) and 



^(^(^A)), 



•(^ ± (x,A)) j± ^(A) 



with(j+,j ) = (l,...,fc, fc + l,...,n). 

Observe that for all x G R and A G Q, y o and Zq satisfy 

Z +(x,A)y -(x,A)=/ fe , 



(15a) 
(15b) 
(15c) 
(15d) 



Zq(x, A)F + (x, A) = I n - k , 
Z£(x, A)F + (x, A) =0(„_ fe )xfe, 

Z Q (x,\)Y ~(x,\) =0kx(n-k)- 

Similarly, the eigenvalues for the first order system are obtained from 

det F (I L 2 + K (X)R(x)) = 0. (16) 

The domain il mentioned in Hypothesis l2.il is essentially the resolvent set of the operator Lq, that is 

fl = C\o- e (L), (17) 

where o~ e (L) = er e (Lo) is the essential spectrum of L given by 



o- e {L ) = {\eC: A = ^c;C': (eK, Q = i l ai and c„ = i n }, 

1=0 

or equivalently 

o- e {L A ) = {A G C : det (i(I n - A(X)) = 0, C G »}• 

Note that A = is an eigenvalue of the operator L embedded in the essential spectrum. 

3 The Birman-Schwinger operator 

Since our problem of computing the pure point spectrum of the operator L is reduced to the computation of 
the determinant of a Fredholm operator given by an identity plus a compact operator. Then, we denote by 
C(%) the set of compact operators in a separable Hilbert space H — L 2 , and by a n (A) the singular values of 
a compact operator A. The Schatten-von Neumann classes of compact operators are defined as 



J P {H)= {AeC{n): 



-|i/p 



with the corresponding norm 



\\M P 



X>m 



E ff M 



< oo > , (1 < p < oo) 



1/p 



and the usual convention for p — oo. We say A is of trace class if A G Ji, and is Hilbert-Schmidt if A G Ji. 
If A G Jp, we denote by det p the p- modified Fredholm determinant [52] ■ Then the determinant of the 
Fredholm operator /% + A is given by 



detf (/« + A) = det p (/« + A). 



Assume that v(x) in © is given by 



*)-± (!)&£,. m 



K ml dx k dx T 



for some fixed m € {0, .. .,n — 2}, <j> satisfying (J3J) and Lq a constant coefficient differential operator given 
in ([3]) . Assume Hypothesis 12.11 then for A G O defined in (|17p , equation ((SJ) is written as 



u(x, A) = - (if(A)™) (a;, A) = - / g(x, £, X)v(t)u(£, A) d£, (19) 

ii 

where g(x,£, A) is given in ([?])• 

3.0.1 The case when m = 

When m = 0, v(x) = (f>(x). Let (f>(x) = |0(x)| 1 / 2 e larg( ^^^, then equation (Tl9|) is equivalent to 

^,A) = -(S(A)V)(x,A), (20) 

where %j)(x,X) = |</)(:r)| 1 / 2 e 4ars< ' c ^ :r * l * l u(£, A), and B(X) is the Birman-Schwinger operator 

B{\) = 4>{x)K{\M 1 l\ (21) 

with integral kernel 



l{x,i,\) = < 



£=1 

n 

x) Y, ai(A)e"«"^-O|0(Or /2 , £ < *• 



(22) 



i=fc+l 



Theorem 3.1. |^[ Theorem 4-1 Chap 4] If f € £ p and h E L p with 2 < p < oo 7 i/ien f(x)h(—iV) is in 
J p and 

\\f{x)h{-iV)\\ Jp <{2-K)-Vv\\f\\ L 4h\\ L *. (23) 

For A G $7, we write S(A) = f(x).h(q) = <f>(x).g(q)\(j)\ 1 / 2 where q = 4-, and 



d:r : 
1 



9(q) = — 

J^ajq 3 -A 

3=0 

with g being the inverse of the Fourier transform associated to the constant coefficient differential operator 
Lq — XIjji. Since <j> G L 1 then the function g(x,£, A) given in (l22t is in L 2 in M. x R. Hence 23(A) is a 
Hilbert-Schmidt operator with norm satisfying f|23|) . Furthermore, by observing that g(q) is the product of 
n Hilbert-Schmidt operators and using the argument in [ISJ Theorem 11.2 Chap IV]. It follows that, £>(A) is 
of trace class for A G O. Hence the determinant of the Fredholm operator I L s + B(X) is given by 

detj (I L 2 + 8(A)) = deti (I L 2 + 23(A)) , (24) 

where deti denote the determinant of a trace class operator. Thus we have 

deti (I L 2 + B(X)) = d(X), (25) 

where d(A) is the Fredholm determinant given by 



d(X) = l + J2d n (X) 



with 



d n (X) = — ■■■ det([g(xi,Xj,X)] 

n - Jr Jr 



i,j=l,...,n. 



dx\ ■ ■ ■ dav. 



Since B(X) is of trace class and induced a continuous kernel g(x, £, A) for all (x, £) £ R 2 and AgU. We have 
for n = 1 

k 

di(A) = tr(S(A)) = ^a,(A) / ^(x) dx = ^ a,(A) / 0(») da. (26) 

=fe+i 



y j aj(A) / 0(a;) dx = \^ Qi(A) / </>(a;) dx. 

»=1 ^ R i=k+l ^ R 



3.0.2 The case when ra/0 

Suppose that the integer m € {1, . . . , n — 2}. Then integration by parts of equation ((20)) gives 

u(s,A) = (-l) m+1 [ g m {x,{,\)<f>(t)u(Z,\) d(, 

Jr 

where 

' k 

E (-i) n («r)' n «<We' irM > e<^. 

As in (|2ip . we rewrite (1281) in the Birman-Schwinger form 

^(x ! A) = (-ir+ 1 ( J B m (A)^)( a ;,0, 
where the Birman-Schwinger operator B rn (X) is given by 

B m {\) = ^{x)K m {\)\<p\ l '\ 
with K m (X) the resolvent operator associated to the kernel function given in (J28I) . 



(27) 



(28) 



Corollary 1. For m G {1, . . . ,n — 2}, t/ie operators B m (X) with Bq(X) := £>(A) given in (|21[) are o/ irece 
cZass. 

Proof. Write -B m (A) = 0(x)g m ((7)|0| 1 ' 2 with go := g. From the continuity property of the Green's function 
associated to the operator Lq — A/^2, that is for all m <E {1, . . . ,n — 2}, 



£ a,(A)K-) m = ^a i (A)K+r. 

i=fc+l i=l 

We have 

g m (q) = c m (X)g(q), 

where c m (A) is a non-zero constant. Then it follows that B m (X) is of trace class. 
Theorem 3.2. Suppose that m = n — 1. then the operator J5„_i(A) with kernel given by 

k 



x 



g n -i(x,£,X) = < 



Y / (-^) n - 1 (4) n - 1 ^(X)e K t(^\ ( j > (0\ 1/2 , x<£, ^e 



n 

(*) £ (-ir- 1 («r)«-ia i (A) e «r(-O| 0(e) |i/2 ) ?<x 



D 



is o/ trace class. 



Proof. For m = n — 1, we have -B„_i(A) = (j>(x)g n -i(q)\(j>\ 1 ' 2 . The kernel function <7„_i(x,£, A) is in 
L 2 (R 2 ,C), so B„_i(A) is Hilbert-Schmidt. Observe that 



k 



fJ 



ai(\) ^ (Xi(X) 



(0 = EtFVE 



r~f «7" - i( -u^ K 4 ~ *C 

1=1 * ' J=fc+1 * 



A' 



OiWOtf)*- 1 ^ OiCA)^)" -1 






Then it follows that there exists a non-zero constant c„_i(A) such that 

£n-l(C)=Cn-l(A)5(C)-/*(C), 

where 

MC) = c„. l( A) £ ^7+ E ^ . 

i=k+l l s i=k+l » s 

Hence 

Bn-i(A)-c„-i(A)6(A) + B(A), 

where B(X) = 0(x)ft.(<7)|0| 1 / 2 . The operator B(X) induced a smooth kernel in L 2 (R 2 h ), hence it is of trace 
class. Moreover B(X) is of trace class, then it follows that B„_i(A) is of trace class. □ 

From Corollary [TJ and Theorem 13. 2[ we thus have the following result. 

Proposition 3.3. If B(X) is of trace class then B m (X) is also of trace class, for m G {1, . . . ,ji — 1}. 

3.1 The first order system case 

Under Hypothesis 12.11 and A G il, equation (Tl6l) is 

Y(x, A) = (Ko(X)Y) (x, A) = - / G(x, £, \)R(t)Y(£, A) df, 

where G(a;, £, A) is dehned in (fT4")l. Equivalently, equation (|16p is written as 

*(s,A) = -|i?(.T)| 1/2 (k (\)R*) (3, A) 

where ^(a;, A) = \R(x)\ 1 ^ 2 Y(x, A) and i?(ir) = [/^(x)! 1 / 2 with C/ a unitary transformation so that R(x) — 
R(x)\R(x)\ 1 ^ 2 . Let /C(A) denote the Birman-Schwinger matrix operator defined by 

K(X) = \R(x)\^ 2 K (X)R 
with integral kernel given by 

h( f x ,f-\R(x)\ 1 / 2 Y -(x,X)Z+(^X)R(O, *<£, !,(£!; 

For A G fi, /C(A) is a Hilbert-Schmidt operator. 

Proposition 3.4. For X E Q, the Birman-Schwinger matrix operator /C(A) is of trace class. 



Proof. Write the kernel k as follows 

k[x ^' A) -{k 21 (x,^X) k 22 (x^,X)J 



where fcy(a;,£, A) E L (R , C™**™ 3 ') for i, j = 1,2, ri\ = k and n 2 = n — k. Suppose to = 0, then given the 
structure of the matrix- valued function R(x) defined in (1121) . we have 



*(x,A) 



o 



where -0(a;, A) is introduced in equation (|2T)]). Also fci2 = &2i = &22 = for all (x, £) E M 2 and Aetl. Hence 

k(x, f , A) = diag (g(x, £, A), 0, . . . , 0) . 

This yields /C(A) = diag (B(A), 0, . . . , 0). Thus 

tr(/C(A)) =tr(B(A)) 

and 

det F (7 L 2 + /C(A)) = dct F (I L 2 + 23(A)) = d(X). 

However for to ^ 0, using a transformation essentially integration by parts in the original scalar integral 
equation, lead to the case when m = that is, /C(A) = diag (B m (X), 0, . . . , 0). Thus /C(A) is of trace class, 
for m E {0, . . . , n — 2} and also when m = n — 1. □ 



Given Proposition ^. 41 we have 

trc» (fc(a;, a;, A)) dx = / g m (x,x, A) da; = tr(B m (A)). 



On the other hand, using equation (|17l) . the invariance of trace for a cyclic rotation of three matrices, in 
particular tr(Z^ RY^) = tr(Y" ~ ZqR) and observing that trc»(i?(a;)) = 0. We have 



d;E 



tT C «{k(x,x,X)) dx= / tr C " (\R(x)\ 1/2 Y Q + (x,X)Z Q -(x,X)R(x)) 
Jo ^ ' 

- f tr c „ (\R{x)\ 1/2 Y -(x,X)Z+(x,X)R{x)) dx 

J — oo 

tro (\R(x)\ 1 / 2 (I n -Y -(x,X)Z+(x,X))R(x)) dx 
tr c » f|i?(a;)| 1/2 y" ( r(a;,A)Z+(a;,A)-R(a;)) da; 
tr c „ (±y ± (x,A)^(a;,A)i2(x)) dx. 
Hence the trace of the trace class operator /C(A) is explicitly given by 

tr (/C(A)) = f tr c „ (±1^(3!, A)Z T(a;, A)i?(a;)) dx. 



(30) 



Remark 1. Note that although the kernel k(x, £, A) associated to the trace class operator /C(A) is discontinuous 
along the diagonal. A formula for computing its trace in terms of its kernel function is available. When 
trn>(-R(a:)) ^ Va; E K and /C(A) is trace class then equation (|3"0"|) remains valid (cf. 13 ]). 



4 Fredholm determinant and Evans function 

In this section, we show that the Fredholm determinant associated to the first order system or the scalar 
problem can be reduced to a finite dimensional determinant under the assumption of all the previous hy- 
potheses. Let A fc (C") denote the kth exterior power of the vector space C™ where A is the usual wedge 
product. Giving U = Ux A . . . A U k € A fe (C") and V = Vi A . . . A 14 £ A fc (C"), then the inner product of U 
and V on A (C n ) is defined by 

[U,V] k = deb & .((V i >V j )a>) i,j = l,...,k. 

We introduce the Hodge star operator which is an isomorphism between A n_fe (C n ) and A fe (C"). 
By * we denote the Hodge star operator, so for any U £ A fc (C") and V £ A"~ fc (C™) we have 

lU,*V\ k V = U/\V, 

where V = e\ A . . . A e„ G A n (C n ) is a volume form, and [e.i)i=\... n a unitary basis for C n . 

Definition 4.1. We define the matrix transmission coefficient D(X) by 

D(X)= lim (Yr( x ,\),(Z+(x,\))*)c«, fori,j = l,...,k, (31) 

where 

Y l -(x,X)=(Y -(x,X)) l + f H(x^ 1 X)R(OY l -(t^)d^ where I = 1, . . . , k 

J — oo 

is the Jost solution with H(x,£,X) = Y ~(x,X)Zq(£,X) + Yq~(x, X)Zq(£, A). 

In the context of the Evans function, A is an eigenvalue if the subspaces defined by the set of solutions 
decaying at ±oo of equation (TlT)|) are linearly dependent. Thus, equation (J2J) means that since F + and Zq 
are orthogonal for all x £ K, and in particular when x — > oo. Then, if Y~ and Zq are orthogonal for some 
A £ Q, as x goes to +oo, we must have Y~ and Y^~ are linearly dependent. In other words, if Y~ and Y + 
are linearly dependent for some A £ ft, so are Y~ and Yq when x tends to +oo. 

Lemma 4.2. Let 

U~(x,X) =Y 1 -(x,X)A---AY-(x,X) £ A fe (C") 

and 

S+(x, A) = {Z+(x, X))l A • • • A (Z+(x, X))l £ A fc (C"). 

For X £ CI, we have 

det ck D(X)= lim \U~{x, A), S+{x, A)],. 

a:— )-+c>o 

Proof. Using equation (|12p , we have 

(Yr( x , A), (£+(*, A))*> c » = (s tj + f (Z+& A)) fl(0*r(& A ) d ^ • (32) 

\ J-oo / i t j=\,...,k 

Observe that the integrand of equation (|3"12|) is independent of the variable x, so we have 

dct cfc .D(A)= lim (det ck (Yr(x,X),(Z+{x,X))*} c A 
= lim [l7-(*,A),5^(a;,A)] fc . 

x— > + oo 

n 
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Theorem 4.3. For A £ £1, the determinant of the matrix transmission coefficient equals to the Fredholm 
determinant, that is 

Act ck D{\) =d{\). 

Proof. Taking the limit as x — ► +00 in (|32[) . we have 

det c , D(X) = dct c . Uj + J (Z+(£, X)). R(l;)Yr(S, A) d^ where i,j = 1, • • • , k. (33) 

Observe that equation ([33| can be written as 

det cfc D(X) = det cfc (I k + D(A)) , 

where T>{\) : L 2 (R, C ) — > L 2 (M,C fc ) is an integral operator with kernel h(x,y,X) which is independent of 
the variable x, that is h(x,y,X) = Z^{y, X)R(y)Y~(y, A). The finite rank operator V(X) is analytic in A 
hence (cf. [22]) det^k D(X) is analytic as well. Since h(x,y,X) is continuous, we have (cf. [18] ) 

det ct J(A) = l + y a(P , (A)) (34) 

T=\ 

WllCrC 

k 

C r (D(X))= V / ... / det(/i ipi(t (xp,Xg,A))^ dxi---dx r . 

i 1 ,i 2 ,...,i„=i-' R •« 

Expanding (J33J), thereafter substituting the Neumann series for the Jost solution Y" ; ~(x,A) and reordering 
the multiple integrals, we have 

f °° 1 — 

det cfc £>(A) = 1+ / tro (Y Q -(x, X)Z+(x, X)R(x)) dx + V -C r (V(X)), (35) 

where T>(X) is an integral operator with kernel fc(x,£, A) satisfying 

tr c » (fc(x,x,A)) dx = / tr c „ (±Y±(x, X)Z*(x, X)R(x)) dx. (36) 

Note that the minus sign in the equation (1361) is obtained by using an equivalent definition of -D(A), that is 
D(X)= lim (y+(x,A),(Z -(x,A))*)t>, ior i,j = k + 1, . . . ,n - k. 

x — > — 00 J 

Observe that the Fredholm determinant d(X) associated to the trace class operator /C(A) and detc* D(X) 
are both analytic and their zeros coincide [33]. Furthermore equations (|3l)l) and ([3H]) imply that fc(x,£, A) = 
fc(x,£, A) for all (x,£) S M 2 and X E ft. It then follows that the right-hand side of equation (1331) coincide 
with the expansion of deti [Ijji + /C(A)), that is 

dct cfc D(X) = deti (I L 2 + /C(A)) . 

□ 

Note that from the proof of Theorem 14.31 we conclude that the integral operator T)(X) is of trace class with 
trace given in (|36l) . and that equation ([351) is indeed 

det cfc D(X) = deti (h? + ^(A)) . (37) 

Let us define the Evans function E(X) by 

E(X) = det c » (^~(x, A) r+(x, A)) , (38) 

where Y~ and Y" + are the set of solutions for the first order system ([TO]) that decay at — 00 and +00, 
respectively. Given Theorem 14.31 we have the following. 
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Theorem 4.4. For A G f2, the Evans function and the Fredholm determinant satisfy 



d{X) = m, (39) 



where c(X) = detc-i (Y (x, A) F + (x, A)) 



Proof. Under Hypothesis [2J] det{> (V (a:, A) F + (x, A)) ^ for all A e SI and x G K. Using the Jost 
solutions Yf and Y. + associated to the equation (fT0|) . and observing that from equation (fT5|) . we have 

(Y - V)^ (J). Then 

£(A) = det c „ (Ff(x, A) • • • Yj~(x, A) F+ +1 (x, A) • • • F+(x, A)) 

= dct c „ (F -(x,A) F+(x, A)) det c „ x „ (j n + (J°-[^J) {J~Y-(x,X) J+Y+{x,\)j\ , 
where 

/•rfcoo 

J ± Y ± (x,X) = T H(xA,X)R(OY ± (^X)d^ 

with H(x : ^X) given in Definition l4.il Using equation (fT51) , we have 

m , rmdpt /7 1 * \,(£ 00 ZZ{S,X)R®Y-&X)dt -CZ+(tX)R(OY+(tX)dC\ 

(40) 
Observe that the integrands in the right-hand side of equation (|4U)) are independent of the variable x, so 
taking the limit for instance as x goes to oo, we have 

E(X) = c(A)det c * (l k + ( Z+ {x, X)R{x)Y- (x, X)dx\ . (41) 

Hence by Theorem 14.31 equation (|3U)) follows □ 

Since Lq — XI L 2 is a constant coefficient differential operator then we can identify its determinant with c(A), 
defined in flS^). That is, 

det F (L - XI L 2) = det c » (Y -(x, A) K+(x, A)) . (42) 

Given equation (fBTJf and (|42l) . it seems natural to define the determinant of the Fredholm operator [T71 
Chap XI] Lo + v(x) — XIi,2 to be the Evans function 

E(X) := det F (L Q + v(x) - XI L *). (43) 

Thus, the Fredholm determinant of a trace class operator can be defined as a ratio of the determinants of 
elliptic operators which are Fredholm. In other words, we can extend the notion of the Fredholm determinant 
to non-compact operators. That is, 

, , ,-r ,w>\ / xx dct F (L + v(x) - XI L 2) 

dcti (I L 2 + K(X)v(x)) = \ ' L ' . 44 

det F (L - XI L 2) 

Indeed, let det^ denote the zeta-regularised determinant (cf. [35]). The determinant of an invertible elliptic 
operator A with positive order and admitting a ray of minimal growth is given by 

( d 

det c (A) = cxp I -— (a{s)\ s =o 

with 

U(s) = tr(A- s ) = ]T X- s , 
x n e<j d (A) 
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where <Td{A) denote the discret spectrum of A and Re(s) ^$> 0. Observe that if an elliptic operator A is 
Fredholm then we can define 

det F (A) := dct c (A). 

It then follows that [5S1 Lemma 2.1] 

det F (L + v(x) - XI L 2) = det F [(L - XI L a)(I L 3 + K(X)v{x)] 

= det F (L - A/ L 2)deti(/ L 2 + K(X)v(x)). (45) 

Idcntifiying equations (|4"5"j) and ([55)1 . equation (|4"3")l follows. Observe that away from the ray of minimal 
growth, the multiplicative property (f3T)| of the determinant by means of the Evans function still holds. 
Moreover, the Evans function can be defined when ad(L) = {0} (cf. [2]). Next, we are going to show that 
given equation (|4"Tj) , we can define the p-modified Fredholm determinant as in equation (|4"4")l . 

4.1 A general first order system 

Now suppose that the operator /C(A) is not derived from a scalar problem, and its kernel is defined as in ([2U)) 
with R(x) satisfying 

Pile— e L\ 

tro (R(x)) dx = 0. (46) 

From equation (|41l) . it is clear that the connection between the Evans function -E'(A) and det^kxk D(X) holds 
independently whether the operator /C(A) is of trace class or not. Consequently, equations (J3"S")) and (pi?) 
hold as well. We have the equality of the determinant if and only if /C(A) is a trace class operator. Assume 
that JC(\) is Hilbert-Schmidt. Then for A £ 51 such that ||/C(A)||2 < 1, the Plemelj formula for det p and in 
particular for p = 2 is given by 

/ °° f_iy-i \ 

det 2 (I L 2 + /C(A)) = cxp I J2 ^— tT (*'(*)) J ■ 

On the other hand, for A £ 51 such that ||X>(A)||i < 1, we write 

det c * D(X) = deh(I L 2 + 35(A)) = cxp (j^i^lj tr (p'(A)) j . (47) 

We know that if /C(A) £ Ji then by Theorem[01 we have tr (35* (A)) = tr (/C'(A)) for I > 1. However /C(A) 
is Hilbert-Schmidt then its trace might not be defined or might not be equal to J R trc« (k(x,X, A)) da;. In 
other words, we have tr (35* (A) J = tr (JC 1 (X)) for I > 2. Thus from (gTJ), we have 

det £fc D(X) = exp (tr (35(A))) det 2 (/ L 2 + /C(A)), (48) 

where tr (35(A)) is defined in (l36l) . Combining (|48|) and (l4lT) . we have 

dct 2 (/ L2 + K(X)) = ^ exp (-tr (15(A))) . 



Consequently, from (J4TJ) and (1431) we have 

More generally if /C(A) € J/^ for p > 2 we have 

det cfc £>(A) = exp j J^ ^—J tr (^( A )) ) det P i 1 ^ + K(A)) ■ 
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Hence 

det p (I L i + /C(A)) = B& exp ( ^ ^-tr [W(X) ) ) . (49) 

Remark 2. Note that equation (I46[) is an immediate consequence of (|36[) . In other words, to ensure the 
equality of the k and n — k determinant in (j4T"j) . equation ([4"o"[) must be satisfied. The condition (14^)) is 
equivalent to the condition J trc n (A(x))dx = given in [13, Theorem 3.2]. 

In the travelling wave problem, it may happen that a transition to instability occurs (see [51 134] ) due to an 
eigenvalue emerging from the essential spectrum. For that case, an extension of the Evans function across 
the essential spectrum is possible (see [2(51 [T2] ) . However, given ([4*9")) we can also use the p-modified Fredholm 
determinant to locate these values that become eigenvalues. Observe that the subdomain $7 can be extended 
to a domain T in the essential spectrum in which c(A) =/= for all A G C. In particular, if the essential 
spectrum is in the left of the spectral plane, we have 

T = {Rc(A) < : j~J {m - Kj) ^ with Rc(kj) G R}. (50) 

l<i<j<n 

Note that for A G T, (Lq-XI^)^ 1 is not an operator in L 2 (R). However if G £ X (R, e^^da;) for some (3 > 
then we might have (f>(x)(Lo — A/^ 2 ) _1 (S)2/)|0(y)| 1 ^ 2 in L 2 (W 2 ). So its induced operator is Hilbert-Schmidt. 
Consequently det2 (Ijji + /C(A)) is analytic for A G T and its zeros are precisely these values that become 
eigenvalues. 



5 Fredholm determinant for fronts 

In this section, we construct an integral operator associated to the Fredholm determinant for the case of 
a front. From previous sections we know that computing the Fredholm determinant corresponding to the 
scalar or the first order problem is the same. Then we can just focus on the first order system case. 

Hypothesis 5.1. Suppose that lim^^ioo R(x) = R , where R + / R~ are constant-valued matrices , and 
\\R(x) — R Hc^xn 6 L 1 satisfying 

/0 p+oo 

tro (R(x) - R-) dx + / tr C n (R(x) - R+) dx = 0. 
-oo J0 

Moreover, assume that the equation 

det c ^ n (A ± (\)~fiI n ) =0, 

where A (X) — A(X) + R , has k roots n i — ii i (A), i — 1, ...,k with Re(K i ) > and n — k roots 
Tf 1 = ^(X), i = k + 1, . . . , n with Re(r i ) < for all X 6 ft, where ft C C is the region of interest. 

When lim-E^y-i-oo R(x) = R , the construction of the Green's matrix associated to the unperturbed problem 
d/dx — A ± (X) is not obvious as in Section [2J However, we observe that we can construct an equivalent 
constant-coefficient differential operator d/dx — B(X) which retain only important informations about the 
operator d/da; — A ± (X). To this end, let us construct the matrix B(\) such that the functions (Yq)j G 
Ker (d/dx — A ± (X)) and decaying at ±oo, satisfy also 



_d 
dx 



T -(Xo ± )j = B(x)(y Q ± ) J . 



It then follows that B(X) must satisfy 

P- 1 (X)B{X)P(X) =diag( K f 
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where 



/ 1 



P(X) = 



'k+l 



VCT)"" 1 



i ^ 



(r+) B "V 



Thus the Green's matrix corresponding to the operator d/dx — -B(A) is given by 

' -Y -(x,X)Z+(£,X), x<C, x,CeR; 



G(z,£,A) 



where 



r +(x,A)Z -(£,A), £<x, 
A)) j = -(Z ± (x,A)) i i?(A). 



(51) 



{Z±(x 



dx V u V ' v yj v u v- • v,j 

The Green's matrix constructed in f|5 1 [) satisfies equation (|12p. In the same manner, the Birman-Schwinger 
operator corresponding to the new problem 

d 



da; 



Y = (B(X) + Q(x))Y 



is given by 



^.■(A) = K t )(, : )| 1 /a(^_B(A) 



Q(x) 



where 

'r{x)-R-, x < 0; 

R(x)~R+, x>0, 

and Q(a;) = {/^(a;)) 1 / 2 with [/ a unitary transformation so that <3(x) = Q(x)|<5(a;)| 1 / 2 . 

Suppose that the Birman-Schwinger operator associated to the original problem, /C(A) and ICb(X) are in J p . 

Then we have 

det p (I L 2 +K(X)) = c(A) det p (I L z + K B (A)) , (52) 

where c(A) is a non-zero analytic function. If for instance /C(A) £ J\ then by Theorem 14.31 we have 



det c „(r-(.x,A) Y+(x,X)) 



c(A)deti(/ L 2+/C B (A)) 



det C n (Y -(x,X) Y+(x,X)) 
The subregion fl when R~ ^ R + is given by [20J Lemma 2, p. 138] 

n = C\(<r e (L A -)Ua e (L A+ )), 
where L A ± = d/dai — A ± (A) . 

Remark 3. If the Evans function is given by (cf. [1]) 

23(A) = cxp ( - J trA(s, A)ds J detc» (V _ (a:, A) Y+(x, A)) . 

Then using Proposition 2 and an argument found in [3], we can show that 

£(A)=C(A)det c ™(Y-(x,A) F+(a;,A)), 

where C(X) is a non-zero analytic function. Hence the connection of E(X) and the p-modified Fredholm 
determinant follow from the previous section. 

Remark 4. Given the connection of the Evans function and the p-modified Fredholm determinant in the one 
dimensional case, a connection in a higher dimensional problem might be possible. This is either accomplished 
by converting the higher dimensional problem to a one dimensional problem using a suitable transformation 
(see. [15]). for numerical details (cf. 31 ). Thus, under hypotheses described in this paper, the connection 
follows. Or, it might be accomplished by constructing the generalised Evans function (see. [B]). 
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6 Conclusion 

We have investigated the connection between the Evans function and the p-modified Fredholm determinant 
for the travelling wave problem in one dimension. This connection turned out to be natural, therefore 
our next objective is to numerically compute the p-modified Fredholm determinant using two approaches. 
These approaches are, the Nystrom method which involves the integral operator itself (see. [3]) and the 
shooting method which uses the Evans function see for example [23J 1301 |3T] for details. Another approach 
is to compute the right-hand side of equation (|4"8")l which involves computing the determinant of a k x k or 
n — k x n — k matrix. 
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